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ABSTRACT 

Why does the Zel'dovich approximation (ZA) work well to describe gravi- 
tational collapse in the universe? This problem is examined by focusing on its 
dependence on the dimensionality of the collapse. The ZA is known to be exact 
for a one- dimensional collapse. We show that the ZA becomes progressively more 
accurate in the order of three-, two-, and one- dimensional collapses. Furthermore, 
using models for spheroidal collapse, we show that the ZA remains accurate in 
all collapses, which become progressively lower dimensional with the passage of 
time. That is, the ZA is accurate because the essence of the gravitational collapse 
is incorporated in the ZA. 

Subject headings: cosmology: theory — large scale structure of universe 
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1. INTRODUCTION 

For gravitational collapse in the universe, Zel'dovich (1970) proposed an approximation, 
the Zel'dovich approximation (ZA). The ZA is a Lagrangian approximation based on the 
motion of a fluid element. If its initial position is q in comoving coordinates, the Eulerian 
position x(q,t) at time t is 

x(q,t) = q + *(q,t). (1) 
The displacement ^(q, t) satisfies the Lagrangian equation of motion (Peacock 1999, §15.2): 

_ + 2ff _._V.*. (2) 

Here H(t) is the Hubble parameter. The gravitational potential <&(x,t) is obtained from the 
density contrast S(x, t) = det(dxi/ 'dqj)' 1 — 1 and the density parameter Q(t) via the Poisson 
equation 

V> = |ff 2 fi<5. (3) 

From the first-order perturbation to equations (2) and (3), Zel'dovich (1970; see also Shan- 
darin & Zel'dovich 1989; Sahni & Coles 1995) assumed that the displacement * scales with 
the linear growth factor D(t): 

2D 

* = " wHT VA ' (4) 

Here the subscript "in" indicates a value at the initial time t in . The linear growth factor D 
satisfies the evolution equation 

Hence, the displacement * is assumed to satisfy the following equation: 

+ 2H^- = -H 2 ^. (6) 



dt 2 dt 2 

Equation (4) or (6) provides an approximation for the gravitational collapse, until the shell 
crossing where fluid elements from different initial positions reach the same position in Eu- 
lerian coordinates. The approximation applies to any cosmological model. 

The ZA is used for various topics in observational cosmology, e.g., the construction of 
the initial conditions in an A"-body simulation (Doroshkevich et al. 1980; Bertschinger 1998), 
the estimation of the initial conditions of the universe from observations of galaxies in the 
present universe (Dekel 1994), and the analytical study of the pairwise peculiar velocities of 
galaxies (Yoshisato et al. 2003). The ZA is also used to develop the adhesion approximation 
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(Gurbatov et al. 1989), the truncated ZA (Coles et al. 1993), and the post-, post-post-, 
and Pade-ZAs (Bernardeau 1994; Buchert 1994; Munshi et al. 1994; Bouchet et al. 1995; 
Catelan 1995; Matsubara et al. 1998). Thus, the ZA plays an important role. 

Despite its simplicity, the ZA is accurate even in the quasi-nonlinear regime, provided 
that the shell crossing is unimportant as at early times or at large scales. 1 This accuracy has 
been demonstrated in statistical comparison with an iV-body simulation (Coles et al. 1993; 
Kofman et al. 1994) and with a high-order Eulerian perturbation theory (Bernardeau et al. 
1994; Munshi & Starobinsky 1994). 

Why is the ZA so accurate? The usual explanation invokes the fact that the Lagrangian 
description is intrinsically nonlinear in the density field (Catelan 1995; Bernardeau et al. 
2002). The density contrast 5 is obtained from the displacement Even if the displacement 
is small and lies in the linear regime, the corresponding density contrast could be large and 
lie in the nonlinear regime. 

However, the ZA's accuracy depends on the shape of the overdense region. We believe 
that this dependence is the most essential (Yoshisato et al. 1998; Peacock 1999, §15.8). The 
past comparison of the ZA with an exact solution was made only for a spherical collapse 
(Munshi et al. 1994; Sahni & Shandarin 1996). Yoshisato et al. (1998) made the comparison 
for spheroidal collapses, which are more usual situations in the real universe. It was found 
that, for a given value of the density contrast, the ZA is more accurate in a lower dimensional 
collapse. 2 

Zel'dovich (1970) noticed that, in the ZA equation (6), the acceleration term — 
in the true equation of motion (2) has been replaced with the linear term for the displace- 
ment 3H 2 il^/2. This gravitational acceleration of the ZA is identical to the true gravita- 
tional acceleration in a one-dimensional collapse (Buchert 1989), for which the ZA is exact 
(Doroshkevich et al. 1973). 

Here we analytically show that the ZA's gravitational acceleration approaches the true 
gravitational acceleration in the order of three-, two-, then one-dimensional collapses (§2). 
We also numerically show that any spheroidal collapse becomes progressively lower dimen- 
sional with the passage of time. The ZA's gravitational acceleration remains close to the 



1 Yoshisato et al. (2003) proposed a method for determining the scale range where the shell crossing is 
unimportant at a given time. 

2 Throughout this paper, as in Yoshisato et al. (1998), we use the dimension to characterize the shape of 
the collapse, which is assumed to occur in a three-dimensional space. The one-, two-, and three-dimensional 
collapses respectively correspond to planar, cylindrical, and spherical collapses. 
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true gravitational acceleration (§3). Then we discuss why the ZA is accurate (§4). 



2. DIMENSIONAL DEPENDENCE OF EQUATION OF MOTION 

The Lagrangian equations of motion (2) are derived for one-dimensional planar, two- 
dimensional cylindrical, and three-dimensional spherical collapses. We compare these equa- 
tions with the ZA equation (6). 



2.1. One-dimensional Planar Collapse 



When the collapse is planar, the gravitational potential t) is obtained using Green's 
function 3H 2 il\x\/A for the one- dimensional Poisson equation (3): 



O fOO 

$ = -H 2 n 



/oo 
s(x',t)\ 
-oo 



x — x'ldx'. 



(7) 



The gravitational acceleration of the fluid element at position x is 

<9$ 3 



dx 



--H 2 n 

4 



POO 

5(x',t)dx'- / 5(x',t)dx' 

i J X 



(8) 



This gravitational acceleration is proportional to the difference of the integrated density con- 
trast between both sides of the position x. The reason is that, in a one-dimensional collapse, 
the gravity is independent of the distance. We rewrite equation (8) with the displacement 
of the fluid element \I> = x — q. Here q is the initial position. Until the shell crossing, mass 
conservation leads to 



1 + 5(x', t)dx' 
1 + S(x', t)dx' 



— oo 
oo 



dx' . 



(9) 



Here we set the initial density contrast to zero. It follows from equation (9) that the inte- 
grated density contrast is proportional to the displacement 



/x pq px 

S(x',t)dx' = / dx' - dx' = q - x = 
oo J —oo J —oo 

POO POO 

5(x', t)dx' = / dx' - dx' = x - q = 

J q J x 



oo 
oo 



(10) 
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The gravitational acceleration is 

= 3 -H*m. (ii) 

dx 2 v ; 

Thus, in the one-dimensional planar collapse, the true gravitational acceleration is identical 
to the ZA's gravitational acceleration. 



2.2. Two-dimensional Cylindrical Collapse 

When the collapse is cylindrical, we use Green's function 37J 2 S7 In |cc|/47r for the two- 
dimensional Poisson equation (3). The gravitational acceleration of the fluid element at 
position x is 

8<b 3 1 r\ x \ 

-a^r-2 H ^\L Hr ' t)r ' dr '- (12) 

This acceleration is proportional to |cc| -1 times the density contrast integrated over the radius 
\x\, in accordance with Gauss's theorem (Binney & Tremaine 1987, §2.0). The gravity in a 
two-dimensional collapse is proportional to the inverse of the distance. We rewrite equation 
(12) with the initial position of the fluid element q. Until the shell crossing, mass conservation 
leads to 

'"I r\q\ 
[1 + S(r', t)]r'dr' = / r'dr'. (13) 
'o Jo 

Then, as in §2.1, the gravitational acceleration is 

d * fon^ , ■ (14) 



\x\ 


2 




\q\ 


2 




x\ 





d\x\ 4" 

Since the displacement ^ is \x\ — \q\, this true gravitational acceleration differs from the 
ZA's gravitational acceleration 3H 2 il^/2 by a factor (|cc| + |q|)/2|x|. 



2.3. Three-dimensional Spherical Collapse 

When the collapse is spherical, we use Green's function — 3H 2 il/8irr for the three- 
dimensional Poisson equation (3). The gravitational acceleration of the fluid element at 
position x is 

f)ti> s i r\ x \ 

-70 = -^ n n2/ W^dr'. (15) 

o\x\ 2 \x\ 2 J 

This acceleration is proportional to |cc| -2 times the density contrast integrated over the radius 
\x\, in accordance with Gauss's theorem (Binney & Tremaine 1987, §2.0). The gravity in a 
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three-dimensional collapse is proportional to the inverse of the distance squared. We rewrite 
equation (15) with the initial position of the fluid element q. Until the shell crossing, mass 
conservation leads to 

r\x\ r \q\ 

/ [l + 5(r',t)]r' 2 dr' = / r' 2 dr'. (16) 
Jo Jo 

Then, as in §2.1 and 2.2, the gravitational acceleration is 



d\x\ 2 |a?| 

Since the displacement ^ is \x\ — \q\, this true gravitational acceleration differs from the 
ZA's gravitational acceleration 3H 2 £l^/2 by a factor (\x\ 2 + \x\\q\ + |q| 2 )/3|a;| 2 . 



2.4. Dimensional Dependence 

Summarizing the equations of motion for one-dimensional planar, two-dimensional cylin- 
drical, and three-dimensional spherical collapses (eqs. [11], [14], and [17]), we have 

dt 2 dt 2 Ix^- 1 n y 1 

with the displacement \I/ = \x\ — \q\ and the dimension n — 1, 2, or 3. The factor l/|a3|™ _1 
reflects the dependence of the gravity on the distance |cc|. The factor (|cc| n — \q\ n )/n reflects 
the change of the integrated density contrast in the course of the gravitational collapse. 

We compare equation (18) with the ZA equation (6), in order to study the ZA's accu- 
racy. When the displacement is small, these two equations are almost identical. The ZA is 
accurate. Even when the displacement is large, in a lower dimensional collapse, the gravity 
l/|a?|™ _1 is more independent of the distance \x\, and the change of the integrated density 
contrast (|a;| n — \q\ n )/n is closer to the displacement \x\ — \q\. The ZA becomes progres- 
sively more accurate in the order of three-, two-, and one-dimensional collapses. For the 
one-dimensional collapse, the ZA equation (6) is exact. These results are independent of the 
functional forms of H(t) and Q(t) and thus apply to any cosmological model. 



3. SPHEROIDAL COLLAPSE 



During a real gravitational collapse, the overdense region changes its shape. The effect 
of this change on the ZA's accuracy, especially in the gravitational acceleration, is studied 
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by using the gravitational collapse of a homogeneous spheroid as an idealized model (see also 
Yoshisato et al. 1998). 

When pressure is negligible, the gravitational collapse increases the ellipticity of the 
spheroid (Lin et al. 1965). While an initially oblate spheroid forms a plane-like structure, 
an initially prolate spheroid forms a line-like structure. This is because the gradient of the 
gravitational potential generated by the spheroid is larger along the minor axis of the spheroid 
than along the major axis. The gravitational acceleration is stronger and the gravitational 
collapse is more significant along the minor axis (Binney 1977). Since the collapse becomes 
progressively lower dimensional with the passage of time, the ZA is expected to remain 
accurate. 



3.1. Equation of Motion 



The density contrast in the presence of a homogeneous ellipsoid is determined by the 
half-lengths of its principal axes aiit): 



S(x,t) = 5 c (t)S 



x 2 



X\ 



a\(t) a 2 2 (t) a$(t) 



(19) 



Here Q(x) is a step function: Q(x) = 1 for x > and for x < 0. The density contrast 
within the ellipsoid S e (t) is 



*,(*) = (! + fin) ";^ 2 ;^'" -i- 

a 1 {t)a 2 {t)a 3 {t) 



(20) 



We adopt spheroidal symmetry, a 1 (t) = a 2 (t). Then the Poisson equation (3) yields the 
gravitational potential Q(x,t) as 



i=l 



with 



A 1 (t) = A 2 (t) = |[1 + h(t)}, A 3 (t) = |[1 - 2h{t)\. 



For an oblate case (aii = a 2 > a 3 ), we have 

h(t) = 

with the ellipticity 

e(t) = I 



3v / T^? 

2 1F~ 



•sin 1 e 



3-e 2 
2e 2 ' 



1/2 



(21) 



(22) 



(23) 



(24) 
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For a prolate case {a.\ = a<i < 0:3), we have 3 



Ht) = 7 



31 -e 2 



4 e 3 



In 



l-e 

TTi 



+ 



3 -2e 2 

2e 2 : 



(25) 



with the ellipticity 



e(t) 



_ ai(t) : 
« 3 (t) 2 



2H 1/2 



(26) 



These equations are derived in Binney & Tremaine (1987, §2.3). We consequently obtain 
the equation of motion (2) for the half-lengths oti(t) 



--H 2 Q5 e A iai . 



(27) 



The evolution of the half-lengths cti(t) describes the evolution of the homogeneous spheroid. 



3.2. ZA 

Equations (4) and (21) yield the equation for the evolution of the half-lengths cti(t) in 
the ZA: ] 

ol% - <Xi,m = --D5 cMl A iMl a itin . (28) 

The ZA equation of motion (6) is 

^1 + 2H^ = -H 2 n(a t - a, in ). (29) 
dt 2 dt 2 v ' ; v ; 

The factor is larger along the minor axis of the spheroid than that along the major axis. 
Hence, also in the ZA, the gravitational collapse is more significant along the minor axis. 
The ellipticity e(t) or e(t) accordingly increases. 



3.3. ZA versus True Solution 



Equation (27) is numerically solved to obtain the evolution of the homogeneous spheroid. 
For simplicity, we adopt the Einstein-de Sitter model (Q — 1 and H = 2 /St), but the results 
could qualitatively apply to any cosmological model. The initial axis ratio ai^/a^in or 



3 The eqs. (33) and (36) in Yoshisato et al. (1998) that correspond to our equations (20) and (25) have 
typographical errors. 
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tt3,in/tti,in is 2, 4, or 8. The initial density contrast is 5 Ci i n = D in . The initial velocity is 
determined using the ZA (eq. [28]). We also obtain the evolution in the ZA. 

Table 1 shows the scale factor at the moment of complete collapse a co \. It is evident 
that, in a lower dimensional collapse, the ZA provides a more accurate result. 

Figures 1 and 2 show the results for oblate collapses (a\ = a<i > a^) and for prolate 
collapses (ol\ = a<i < 0:3), respectively, as a function of the normalized scale factor a/a co \ 
(solid lines). For reference, the solutions for the one-dimensional planar, two-dimensional 
cylindrical, and three-dimensional spherical collapses are also shown (dotted lines). 

First, we discuss the evolution of the half-lengths (Figs, la and 2a). The collapse is 
rapid along the minor axis, while it is slow along the major axis. In the oblate cases (Fig. 
la), the collapse asymptotically becomes one-dimensional. The collapse along the minor 
axis is intermediate between the one- and three-dimensional collapses. If the initial axis 
ratio a^in/atyin is large, the collapse is close to the one-dimensional collapse. In the prolate 
cases (Fig. 2a), the collapse asymptotically becomes two-dimensional. The collapse along 
the minor axis is intermediate between the two- and three-dimensional collapses. If the 
initial axis ratio a^m/tti.in is large, the collapse is close to the two-dimensional collapse. 

Second, for the gravitational acceleration along the minor axis, we compare the ZA 
with the true numerical solution (Figs, 16 and 26). In the oblate cases (Fig. 16), the ZA's 
accuracy is intermediate between those in the one- and three-dimensional collapses. If the 
initial axis ratio is large, the ZA's accuracy is close to that in the one-dimensional collapse 
in which the ZA is identical to the true solution. This is because, as shown in Figure la, 
the collapse is close to the one-dimensional collapse. In the prolate cases (Fig. 26), the ZA's 
accuracy is intermediate between those in the two- and three-dimensional collapses. If the 
initial axis ratio is large, the ZA's accuracy is close to that in the two-dimensional collapse. 
This is because, as shown in Figure 2a, the collapse is close to the two-dimensional collapse. 
The relative accuracy is the same if the ZA is compared with the true numerical solution 
along the density contrast S e in Figure 3. 

Third, for the density contrast, we compare the ZA with the true numerical solution 
(Figs, lc and 2c). In the oblate cases (Fig. lc), the ZA's accuracy is intermediate between 
those in the one- and three-dimensional collapses. If the initial axis ratio is large, the ZA's 
accuracy is close to that in the one-dimensional collapse. In the prolate cases (Fig. 2c), the 
ZA's accuracy is intermediate between those in the two- and three-dimensional collapses. 
If the initial axis ratio is large, the ZA's accuracy is close to that in the two-dimensional 
collapse. The ZA is less accurate in the density contrast than in the gravitational acceleration 
at the late stage of the evolution. This is because complete collapse occurs earlier in the true 
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numerical solution than in the ZA (Table 1). 

4. DISCUSSION 

The ZA is the linear approximation to replace — V x $ with 3H 2 n^/2 for the gravita- 
tional acceleration in the Lagrangian equation of motion (eqs. [2] and [6]). This approxima- 
tion becomes progressively more accurate in the order of three-, two-, and one-dimensional 
collapses (eq. [18]). The reason is that, in a lower dimensional collapse, the gravity 1 / 1 a? | ri 1 
is more independent of the distance \x\, and the change of the integrated density contrast 
(|cc| n — \q\ n )/n is closer to the displacement ^ = \x\ — \q\. For the one-dimensional collapse, 
the ZA's gravitational acceleration is exact. 

The spheroidal collapse preferentially proceeds along the minor axis, i.e., the direction 
of the strongest gravitational acceleration (Figs. 1 and 2). While an oblate spheroid asymp- 
totically undergoes a one-dimensional collapse, a prolate spheroid asymptotically undergoes 
a two-dimensional collapse. The dimension of the collapse decreases in the course of the 
collapse. This is well reproduced in the ZA. The ZA's gravitational acceleration remains 
close to the true gravitational acceleration. 

During structure formation in the real universe, an overdense region preferentially col- 
lapses along the direction of the strongest gravitational acceleration. The collapse asymp- 
totically becomes one-dimensional as in the case of oblate spheroidal collapse. 4 

Therefore, the ZA is accurate for the following reasons. First, the ZA is more accurate 
in a lower dimensional collapse (see also Yoshisato et al. 1998). Second, any gravitational 
collapse preferentially proceeds along the direction of the strongest gravitational acceleration, 
and it becomes progressively lower dimensional with the passage of time. That is, the ZA is 
accurate because the essence of the gravitational collapse is incorporated in the ZA. 

This conclusion applies to other approximations based on the ZA, e.g., the truncated ZA 
(Coles et al. 1993). Such an approximation has inherited the ZA's dependence of accuracy 
on the dimensionality of the gravitational collapse (see also Yoshisato et al. 1998). 

Our study is based on Buchert (1989), in which the ZA's gravitational acceleration 
was found to be exact for a one-dimensional collapse. This finding was for local motions. 



4 Since, at a later stage, the overdense region could collapse along the other directions, the final structure 
is not necessarily plane-like but could be point- or line-like (Arnold ct al. 1982; see also Sahni & Coles 1995; 
Kurokawa et al. 2001). The ZA is not accurate, because of the shell crossing, after the complete collapse 
along the direction of the strongest gravitational acceleration. 
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We have shown that this is also true for global motions. We have also shown that the 
ZA's gravitational acceleration is closer to the true gravitational acceleration in a lower 
dimensional collapse. 

The usual explanation for the ZA's accuracy, especially in the density, is that the den- 
sity contrast is a nonlinear function of the displacement (Catelan 1995; Bernardeau et al. 
2002). Even if the displacement is small and lies in the linear regime, the corresponding 
density contrast could be large and lie in the nonlinear regime. For a given change of the 
displacement, the corresponding change of the density contrast is larger in a higher dimen- 
sional collapse (eq. [18]). The ZA should be more accurate in a higher dimensional collapse, 
according to this usual explanation. However, in practice (Figs. 1-3; see also Yoshisato et 
al. 1998), the ZA is more accurate in a lower dimensional collapse. 

The ZA has long been known to be exact for a one-dimensional collapse (Doroshkevich 
et al. 1973). However, with this property alone, it is impossible to explain the ZA's accuracy. 
This is because, in general, gravitational collapse is not exactly one-dimensional. Our study 
has clarified that, even if the gravitational collapse is not exactly one-dimensional, the ZA 
is sufficiently accurate. Thus we have successfully explained why the ZA is accurate. 

The authors are grateful to the referee for careful reading of the manuscript and helpful 
comments. 
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Table 1. Scale Factor a co i at the Moment of Complete Collapse 



Initial Shape 


True 


ZA 


three-dimensional spherical 


1.686 


3 


two-dimensional cylindrical 


1.466 


2 


one-dimensional planar 


1 


1 


oblate 0:1/0:3 = 2 


1.412 


1.897 


Oii/013 = 4 


1.216 


1.421 


Oil/ 0(3 = 8 


1.110 


1.204 


prolate 0:3/01 = 2 


1.577 


2.420 


03/01 = 4 


1.510 


2.163 


03/01 = 8 


1.480 


2.059 
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Fig. 1. — Oblate spheroidal collapse, (a) Half-lengths a,\ and 0:3 in the true solution, (b) 
The ZA's gravitational acceleration along the minor axis 3H 2 Q(as — a3 ) - m )/2 in the ZA. (c) 
Density contrast 5 C in the ZA. The half-lengths are normalized by their initial values. The 
ZA's gravitational acceleration and density contrast in the ZA are respectively normalized by 
the true gravitational acceleration — 3H 2 Q5 e A 3 a 3 /4: and density contrast in the true solution. 
The abscissa is the scale factor a normalized by its value at the moment of complete collapse 
a co \ in the true solution. The dotted lines denote the one- and three-dimensional collapses. 
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a C0 |(true) 



Fig. 2. — Same as Fig. 1 but for prolate spheroidal collapse. The ZA's gravitational acceler- 
ation is 3H 2 £l(ai —a 1: - m )/2 and the true gravitational acceleration is — 3H 2 Q5 e Aicti/4: along 
the minor axis. The dotted lines denote the one-, two-, and three-dimensional collapses. 
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Fig. 3. — The ZA's gravitational acceleration along the minor axis normalized by the true 
gravitational acceleration, (a) Oblate spheroidal collapse, (b) Prolate spheroidal collapse. 
These are the same as Fig. lb and 2b, respectively, but the abscissa is the density contrast 
S e in the true solution. 



